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The impact of the spin-flip terms on the (static and dynamic) charge and spin correlations in 
the Luttinger-liquid ground state of the one-dimensional (ID) t-J model is assessed by comparison 
with the same quantities in the ID t-Jz model, where spin-flip terms are absent. We employ the 
recursion method combined with a weak-coupling or a strong-coupling continued-fraction analysis. 
At Jz/t = 0^ we use the PfafBan representation of dynamic spin correlations. The changing nature 
of the dynamically relevant charge and spin excitations on approach of the transition to phase 
separation is investigated in detail. At the transition point, the t-Jz ground state has zero (static) 
charge correlations and very short-ranged (static) spin correlations, whereas the t-J ground state 
is critical. The t-Jz charge excitations (but not the spin excitations) at the transition have a 
single-mode nature, whereas charge and spin excitations have a complicated structure in the t-J 
model. A major transformation of the t-J spin excitations takes place between two distinct regimes 
within the Luttinger-liquid phase, while the t-Jz spin excitations are found to change much more 
gradually. In the t-Jz model, phase separation is accompanied by Neel long-range order, caused by 
the condensation of electron clusters with an already existing alternating up-down spin configuration 
(topological long-range order). In the t-J model, by contrast, the spin-flip processes in the exchange 
coupling are responsible for continued strong spin fiuctuations (dominated by 2-spinon excitations) 
in the phase-separated state. 



*Present address: C.W. Costello & Associates, Providence RI 02903 



1 



I. INTRODUCTION 

At the heart of many phenomena in condensed mat- 
ter physics is the interplay between the charge and spin 
degrees of freedom of interacting electrons. The impact 
of the magnetic ordering and fluctuations on the charge 
correlations or the effect of the phase separation on the 
spin correlations, for example, are important issues in 
the study of strongly correlated electron systems. One 
of the simplest scenarios in which these questions can 
be formulated transparently and investigated systemat- 
ically comprises two successive approximations of the 
Hubbard model with very strong on-site repulsioa. They 
are known under the names t-J and t-J^ modelsE- _ 

Here we consider a one-dimensional (ID) lattice. In 
both models the assumption is that the Hubbard on-site 
repulsion is so strong that double occupancy of electrons 
on any site of the lattice may as well be prohibited com- 
pletely. This constraint is formally incorporated into the 
two models by dressing the fermion operators of the stan- 
dard hopping term with projection operators: 



Ht = -t 



:,.} (1-1) 



with ci^a = c;,o.(l-n;,_c,), ni = ni^^+ni^i, ni^„ = 
In the t-J model the Hubbard interaction is further taken 
into account by an isotropic antifcrromagnctic exchange 
coupling between electrons on nearest-neighbor sites: 



Ht-J — Ht + J ^ I Si ■ S; + i 



1 



-nm+i 



(1.2) 



with Si = |(ni,T - "-/a); 



ci^l, and Si 



In the t-Jz model the isotropic exchange interaction is 
replaced by an Ising coupling: 



(1.3) 



The absence of spin-flip terms in Ht.j^ introduces ad- 
ditional invariants (not present in Ht-j) for the spin con- 
figurations of eigenstates and thus alters the relationship 
between charge and spin correlations considerably. All 
results presented here will be for one-quarter-filled bands 
(A^e — N/2 electrons on a lattice of A^ sites). 

For weak exchange coupling, both models have a 
Luttinger liquid ground state. For stronger coupling, 
electron-hole phase separation sets in. Phase separation 
is primarily a transition of the charge degrees of freedom. 
Here it is driven by an interaction of the spin degrees of 
freedom, and it is accompanied by a magnetic transition. 
The degree of spin ordering in the phase-separated state 
depends on the presence {t-J) or absence (t-Jz) of spin- 
flip terms in the interaction. 

Detailed information on the charge and spin fluctu- 
ations in Ht-j and Ht-j^ is contained in the dynamic 



charge structure factor Snnil, ^) and in the dynamic spin 
structure factor Szzil,^), i.e. in the quantity 



+ 00 



SAAiq,(^) 



dte"^\Ag{t)A. 



(1.4) 



where Aq stands for the fluctuation operators 

= AT-^^e-^-^'n,, S'^ ^ N'^ J^e-"^^ S[ . (1.5) 

The degree of spin and charge ordering in the ground 
state is also reflected in the equal-time charge corre- 
lation function {nini+m) and spin correlation function 
(5f S'f_^„j) and in their Fourier transforms, the structure 
factors Snn{q) = {nqU-q) and Szz{q) = (S^S^q). 

In the following we investigate the T = charge 
and spin fluctuations of the two models Ht-j and Ht-j^ 
in three different regimes with the calculational tools 
adapted to the situation: the limit of zero exchange cou- 
pling (Sec. H), the Luttinger liquid state (Sec. [II), and 
the phase-separated state (Sec. |lV[) . 



II. FREE LATTICE FERMIONS 



A. Charge correlations and dynamics 



The tight-binding Hamiltonian (Fl) has a highly 
spin-degenerate ground state. The charge correlations 
are independent of the spin configurations and, there- 
fore, equivalent to those of a system of spinless lattice 
fermions, 



(2.1) 



This Hamiltonian has been well studied in the context of 
the ID 5 = 1/2 XX model, 



which, for J± r^^^^ becomes (2.1) via Jordan- Wigner 
transformation.E£lEJ The equal-time charge correlation 
function of Ht (or H^) exhibits power-law decay, 



cos(7rTO) 



1 



{nini+„i) - {ni){ni+„i) = 
and the charge structure factor has the form 

A^ . 



(2.3) 



Snniq) 



(2.4) 
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The dynamic charge structure factor, whic h is equivalent 
to the zz dynamic spin structure factor of (2.2) reads (for 

2Q(uj - 2t sin q)e{At sin(g/2) - to) 



16i2 sm^{q/2) 



(2.5) 
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B. Spin correlations 

The charge-spin decouphng as is manifest in the prod- 
uct nature of the ground-state wave functions of Ht-j^ at 
Jz/t = 0+ and Ht_j at J/t = 0+ was shown to-lead to 
a factorization in the spin correlation functionfllld'til We 
can write 



m+l 



J=2 



(2.6) 



where C(m) = {SfSij^^LS is the correlation function 
in the ground state of a system of Nf, localized spins 
with antiferromagnetic Heisenberg {t-J) or Ising (t-Jz) 
coupling, and 



P{m,j) EE {nini+,riSj,N^), A^m = ^ ' 



i=i 



is the probability of finding j electrons on sites 1,1 + 
1, . . . , l+m with no holes at the end points of the interval. 
This expression can be brought into the form 



l+ml 



-E- 



S{k) 



X - 2A„-i(fc) + A„-2(fc)], (2.7) 



Sik) 



We 



(2.8) 



Dm{k) = ( exp 



-ik ni 



1 is the 



where S{k) for k = (27r /iVe)n, n = 0, . . 
static structure factor for the localized spins, and the 
Dm{k) are many-fermion expectation values, which arc 
expressible as determinants of dimension m + 



D^{k) 



5t 



(1-t-e-*''-) sin[7r(i - j)/2] 



sin[7r(i - j)/2Ne 



i.j—Q,....m 



lnHt.j^ wehaveC(m) = i(-l)™,i.e. Sik)^iNj'i)Sk,^, 
reflecting the (invariant) alternating up-down sequence 
of successive electron spins. Expression (2.7) can then 
be evaluated in closed form: 



l+2nl 



(~i)"n o2 

2^2 li^*' 

i=l 



(2.9a) 



with 



-2((^r^r+2n> + (^f^r+2«+2» (2.%) 



^*=-n(i-i-i 



The Icading-jtcrms of the long-distance asymptotic ex- 
pansion arellj 



42 1 



4^2 
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1 cos sm 
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(2.10) 



with A = 2^/12 exp[3C'(-l)] = 0.64500 . . .. The structure 
of Dm (tt) is very similar.to that of the xx spin correlation 
function of Hxx .EjuEZI Its leading asymptotic term has 
the form (5f 5f+„)xx ~ (AV2V2)m-i/2. 

In Ht-j the spin-flip terms weaken the spin correlations 
at J jt = O"*". The function S(k) in ( p7^ ) is determined via 
( ^.8[ ) by the spin correlation function of the ID s — 1/2 
Heisenberg antiferromagpet {XXX model). Its leading 
asymptotic term readstSi Cirn) ~ r(— l)'"m~^(lnTO)^/^ 
with amplitude F ~ 0.125(15) as estimated from finite- 
chain data.tj The leading asymptot ic te rm of the t-J pfjin 
correlation function inferred from (|2.7|) has the forrrO 



{Si Si 



l+m) t-J 



TA^V2cos{TTm/2) 



(Inm 



,1/2 



,3/2 



(2.11) 



The t-J and t-J^ spin structure factors Szz{<l) inferred 
from the results presented here will be presented and dis- 
cussed in Sec. III.E. 

For an intuitive understanding of the q = t: charge den- 
sity wave in the ground state at Jz/t = 0+ and J/t = 0+, 
we note that the hopping term opposes electron cluster- 
ing. In the absence of the exchange term, which favors 
clustering of electrons with opposite spin, the hopping 
effectively causes an elect ron repulsion. This is reflected 
in the power-law decay (2.3) of the charge correlation 
function, specifically in the term which oscillates with a 
period equal to twice the lattice constant {q = Akp = t^)- 
In this state, an electron is more likely to have a hole 
next to it than another electron. 

How does this affect the spin correlations? Recall that 
the ground state of Ht-j^ at Jz/t — 0+ is characterized 
by an (invariant) alternating spin sequence. In a perfect 
electron cluster this sequence would amount to saturated 
Neel ordering {q = it), but here it is destroyed by a dis- 
tribution of holes. Spin long-range order exists only in a 
topological sense. However, some amount of actual spin 
ordering survives by virtue of the effective electron repul- 
sion in the form of the algebraically decaying term (2.1C) 
in the spin correlation function with a wavelength equal 
to four times the lattice constant {q = 2kp = 7r/2). 

A similar argument obtains for the t-J model. Since 
its ground state at J/t = 0+ contains all spin sequences 
with S^ = 0, not just the alternating ones, the resulting 
g = 7r/2 oscillations (2.11) in the spin corrd^tion function 
decay more rapidly than in the t-Jz case.E2l 
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C. Spin dynamics 



III. LUTTINGER LIQUID STATE 



Expression (2^) cannot be generalized straightfor- 
wardly for the calculation of dynamic spin correlations, 
the principal reason being that the number of electrons 
between any two lattice sites is not invariant under time 
evolution. However, in the t-Jz case we can determine 
the function {Sf{t)Sf_^^) on a slight detour. We use open 
boundary conditions and write 



1 ' 



(2.12) 



where ctl = ±1 denotes the spin direction of the left- 
most particle in the chain, which is an invariant under 
time evolution. The time-dependent two-spin correlation 
function of the open-ended t-Jz chain is then related to 
the following many-fermion correlation function: 



l+m 



, 

(5f w^f+j = -Ani{t)\{{-ir'^'^ 



c\ {t)ci {t)A^ it)B, (t) A2 {t)B2 (i) • • • Ai {t)Bi {t) 

X A1B1A2B2 ■ ■ ■ Ai+jn Bi^jncj^jn^l+r, 



with A. 



+ ci, Bi 



- ci. In order to extract the 
bulk behavior of {Sf {t)Sf_^„-l) from this expression, we 
must choose both sites I and / + m sufficiently far from 
the boundaries. 

The nurnerical evaluation of this function via Pfaf- 
fians showsclrO that the leading long-time asymptotic 
term describes uniform power-law decay, {Siit)S^^2,i) ^ 
t~^^'^, for even distances and (more rapid) oscillatory 



-2it 



i+2n+l/ 



> 1, 



power law decay, {S[{t)S[^ 
for odd distances. Moreover, we have found compelling 
numerical evidence that the relation (2.9b) can be gen- 
eralized to time-dependent correlation functions in the 
bulk limit / — *■ 00. 

From our data in conjunction with the long-distance 
asymptotic result (2.10) we predict that the leading term 
for large distances and long times has the forrrO 



{^i W^i+m)t-J, ^ 4 (m2- 4*2)1/4 



(2.13) 



The corresponding asjjmptotic result in the XX model 
is well estabhshedMEa 



{snt)si 



1 



A^V2 



XX 



4 (m2 - J2 t2)l/4 



(2.14) 



The asymptotic behavior ( 2.13| ) of the dynamic spin cor- 
relation function implies that the dynamic spin structure 
factor has a divergent infrared singularity at q — n/2: 
S'^2(7r/2,ti;)t_j, ^ Lu~^/^. Further evidence for this singu- 
larity and for a corresponding singularity in Szziq,uj)t-j 
will be presented in Sec. III.F. 



Turning on the exchange interaction in Ht-j and Ht-j, , 
which is attractive for electrons with unlike spins and zero 
otherwise, alters the charge and spin correlations in the 
ground state gradually over the range of stability of the 
Luttinger liquid state. 

In the t-Jz model, where successive electrons on the 
lattice have opposite spins, the exchange coupling coun- 
teracts the effectively repulsive force of the hopping term 
and thus gradually weakens the enhanced q = n charge 
and q — tt/2 spin correlations. We shall see that the re- 
pulsive and attractive forces reach a perfect balance at 
Jz/t = 4~. Here the distribution of electrons (or holes) 
is completely random. All charge pair correlations van- 
ish identically and all spin pair correlations too, except 
those between nearest-neighbor sites. This state marks 
the boundary of the Luttinger liquid phase. At Jz/t > 4 
the attractive nature of the resulting force between elec- 
trons produces new but different charge and spin corre- 
lations in the form of charge long-range order at q = 0^ 
(phase separation) and spin long-range order at g = tt 
(antiferromagnetism) . 

In the t-J model the disordering and reordering ten- 
dencies are similar, but the exchange interaction with 
spin-flip processes included is no longer uniformly at- 
tractive. At no point in parameter space do the attrac- 
tive and repulsive forces cancel each other and produce 
a random distribution of electrons. A sort of balance be- 
tween these forces-jexists at J/t = 2, which is reflected 
in the observationcl that the ground-state is particularly 
well represented by a Gutzwiller wave function at this 
coupling strength. Charge and spin correlations exhibit 
power-law decay at the endpoint, J/t ~ 3.2, of the Lut- 
tinger liquid phase. Here the attractive forces start to 
prevail on account of sufficiently strong antiferromag- 
netic short-range correlations and lead to phase separa- 
tion, but the spin correlations continue to decay to zero 
asymptotically at large distances. 

One characteristic signature of a Luttinger liquid is 
the occurrence of infrared singularities with interaction- 
dependent exponents in dynamic structure factors. In 
the following we present direct evidence for interaction- 
dependent infrared singularities in the dynamic charge 
and spin structure factors of Ht-j^ and Ht-j. We employ 
the recursion methocO in combination with techniques 
of continued-fraction analysis pxpntly developed in the 
context of magnetic insulators .EjI^j 

The recursion algorithm in the present context is 
based on an orthogonal expansion of the wave function 
|*^(i)) = Aq{-t)\(l)) with Aq as defined in (|]|). It 
produces (after some intermediate steps) a sequence of 
continued- fraction coefficients A^(g), A^((7), . . . for the 
relaxation function, 



4 



Co (9,^) = 



z + 



Af(g) 
z + . . . 



which is the Laplace transform of the symmetrized cor- 
relation function The T = dy- 
namic structure factor (|1.4|) is then obtained via 



,)eHlhn3?[c^ 



For some aspects of this study, we benefit from the 
close relationship of the two itinerant electron models 
Ht.j^ and Ht.j with the ID s = 1/2 XXZ model, 

Hxxz = Hxx - J\\ ^ S[S[^i, 



a model for localized electron spins. The equivalence 
of Ht-j^ and Hxxz for Jw^"^ Jzl'^ and J_l = 2t was 
pointed out and used beforeE'Q Depending on the bound- 
ary conditions, it can be formulated as a homomorphism 
between eigenstates belonging to specific invariant sub- 
spaces of the two models. The mapping assigns to any 
up spin and down spin in Hxxz an electron and a hole, 
respectively, in Ht-j^. The spin sequence of the electrons 
in the subspace of interest here is fixed, namely alter- 
natingly up and down. The importance of this mapping 
derives from the fact that the ground statejpmperties of 
Hxxz have been analyzed in great detail.EiTEj 

The r = dynamic charge structure factor Snniq,'-^) 
of Hf.j^ is thus equivalent to the T = dynamic spin 
structure factor Szz{q,^) of Hxxz throughout the Lut- 
tinger liquid phase, and we shall take advantage oL-the. 
results from previous studies of XXZ spin dynamics 0'E£l 
The spin dynamics of Ht-j^ is not related to any known 
dynamical properties of Hxxz- 



A. Charge structure factor 

Certain dominant features of the dynamic charge 
structure factor S'„„(g, uj) are related to known properties 
of the static charge structure factor. Figure |l| displays 
finite- iV data of S'„„(q) for various coupling strengths in 
the Luttinger fiquid phase of (a) Ht-j^ and (b) Ht-j. 

The alignment of the data points on a sloped straight 
line in the free-electron limit represents the exact result 
(2.4), which is common to both models. The persistent 



linear behavior at small q for nonzero coupling reflects an 
asymptotic term of the form Aom~^ in the charge cor- 
relation function (nini^rn), while the progressive weaken- 
ing of the cusp singularity at g = tt reflects an asymptotic 
term of the form ~ AiCos{Trm)/m^p with a coupling- 
dependent charge correlatioft, exponent rjp. For Ht-j^ this 
exponent is exactly knownS 



No exact result exists for the t-J case, but the predic- 
tion is that the charge correlation exponent varies over 
the same range of values,El i.e. between rfp = 2 at J/t ~ 
and On = oo at J/t ~ 3.2. For J/t ^ 1, the data in 
Fig. |l|(b) indicate the presence of a third cusp singu- 
larity in Snn{q), namely at q = 7r/2, which reflects the 
third asymptotic term, ^ A2Coa{nml2)/m^^'^''^^, pre- 
dicted for the t-J charge correlations .l3 No correspond- 
ing singularity is indicated in the data of Fig. |^(a) , nor 
is any corresponding asymptotic term predicted in the 
XXZ spin correlations. 

At the endpoint of the Luttinger liquid phase {Jz/t — 
4), the t-Jz ground-state wave function has the form 



I0o)= E 

l<h<l2<---<lN/2<N 
1 



{n/2 



-1/2 



L, . . . , In/2) 



X -^{inT...)-im...)} , (3.2) 

where \li, ... ,1^/2) specifies the variable charge posi- 
tions. The electrons are distributed completely at ran- 
dom on the lattice, while the sequence of spin orien- 
tations is frozen in a perfect up-down pattern. This 
state is non-degenerate for finite N, and its energy per 
site is A^- independent: Eo/N = —t. For N oo, the 
t-Jz charge correlations disappear completely, {nini+rnj — 
{ni){ni+jn) = as is indicated by the finite-A 

data for Jz/t = 4 in Fig. |l|(a): 5„„(g) - (A/4)5,,o = 
[A/4(A - 1)](1 - (5g,o)- The t-J charge correlations, by 
contrast, seem to persist at J/t 3.2. 



B. Charge dynamics (weak-coupHng regime) 



T]p = 2/[l - (2/7r) arcsin( Jj4i)] 



(3.1) 



Expression ( |2.5| ) for the T — Q dynamic charge struc- 
ture factor Snn (q-, w) of Ht is modified differently under 
the influence of a Jz-type or a J-type exchange interac- 
tion. Within the Luttinger liquid phase we distinguish 
two regimes for the charge dynamics: a weak- coupling 
regime and a strong- coupling regime as identified in the 
context of the recursion methodEj 

In the framework of weak-coupling approaches, the dy- 
namically dominant excitation s pec trum of S'„„((7, w) is 
confined to a continuum as in (|2.5| ) but with modified 
boundaries and a rearranged spectral-weight distribu- 
tion. Moreover, a discrete branch of excitations appears 
outside the continuum. A weak-coupling continued- 
fraction (WCCF) analysis for 5'„„(7r,w) of Ht-j and, in 
disguise, also of Ht-j^ , namely in the form of Szz^t^i ^) for 
Hxxz was reported in Ref. |2^. Without repeating any 
part of that analysis we recall here those results which 
are important in the present context. 

The renormalized bandwidth ujq of the dynamic charge 
structure factor Snni^^, ^) versus the coupling constant as 
obtained from a WCCF analysis is shown in the main plot 
of Fig. H for both the t-Jz model (□) and the t-J model 
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(o). In the XX Z context, wq is the bandwidth of the 2- 
spinon continuum, which is exactly known.E3 Translated 
into t-Jz terms, the expression reads 



ujQ/2t — (n/fi) sin fi, cos/i = — Jz/4i 



(3.3) 



and is represented by the solid line. Comparison with 
our data confirms the reliability of the WCCF analysis. 

Our bandwidth data for the t-J model caii-|be com- 
pared with numerical results of Ogata et alB for the 
charge velocity Vc as derived from the numerical anal- 
ysis of finite chains. The underlying assumption is that 
the relation luq = 2vc, which is exact in Ht_j^, also holds 
for the Ht-j model. The t-J charge-velocity results of 
Rfif. ^ over the entire range of the Luttinger liquid phase 
are shown as full circles connected by a dashed line in 
the inset. The solid line represent s th e exact t-Jz charge 
velocity Vc — with luq from (3.3). 

The dashed line in the main plot is the t-J bandwidth 
prediction inferred from the data of Ref. ^ It is in near 
perfect agreement with the WCCF data (o). The open 
squares in the inset show the WCCF data over a wider 
range of coupling strengths. The renormalized band- 
width loq will shrink to zero at the endpoint of the Lut- 
tinger liquid phase, and the spectral weight will gradually 
be transferred from the shrinking continuum to states of 
a different nature at higher energies. 



J^'^^ ~ 3.2t. Th e short-dashed line represents the exact 
t-Jz result (3.1) thus transcribed for Ht.j. Its deviation 
from the data of Ogata et al. are very small throughout 
the Luttinger liquid phase. 

In Rcf. p9 we carried out a WCCF reconstruction of the 
function S',i„(7r,a;) foe. the t-J model and the t-Jz model 
(alias XX Z model) .EZl The observed spectral- weight dis- 
tributions of both models consisted of a gapless contin- 
uum with a cusp-like infrared singularity (/3p > 0), a 
shrinking bandwidth {uJo/2t < 2), and a lone discrete 
state outside the continuum near its upper boundary. 



D. Charge dynamics (strong-coupling regime) 

What happens to the dynamic charge structure factor 
Snn{q, w) as the exchange interaction is increased beyond 
the weak-coupling regime of the Luttinger liquid phase? 
For the t-Jz case the answer can be infcpjed. from known 
results for the spin dynamics of i/j^xz-l^'E^ The contin- 
uum of charge excitations with sinc-likc boundaries 

T:t sin fl , / \ r. / /r>N 

^L{q)^ \s\nq\, eu{q) ^2eL{q/2}, 

continues to shrink to lower and lower energies, and dis- 
crete branches of excitations 



C. Infrared exponent 

In the Luttinger liquid phase, the dynamic charge 
structure factor has an infrared singularity with an ex- 
ponent related to the charge correlation exponent: 



(3.4) 



The WCCF analysis yields specific predictions for (3p in 
both models. Our results plotted versus coupling con- 
stant are shown in the inset to Fig. js] for Ht.j^ (□) and 
Ht-j (o). The solid line represents the exact t-Jz result 
inferred from (3.1). 

We observe that the WCCF prediction for the infrared 
exponent (□) rises somewhat more slowly from zero with 
increasing coupling than the exact result. The solid line 
in the main plot depicts the inverse square of the exact 
t-Jz correlation exponent (3.1) over the entire range of 
the Luttinger liquid phase. The open squares represent 
the WCCF data for 2 -\- (3p = rip extended to stronger 
coupling. For Ht-j the correlation exponent is not exactly 
known. The solid circles interpolated by the dasied line 
represent the prediction for r]p of Ogata et alJa based 
on a finite-size analysis. The dashed line in the inset is 
inferred from the same data. It agrees reasonably well 
with the WCCF data for /3p (o). 

The solid and long-dashed curves in the main plot sug- 
gest the intriguing possibility that the exponents rjp of the 
two models have the same dependence on the scaled cou- 
pling constants Jzjj'f^ with j'f^ = At and J/J^'^'^ with 



27rtsin^ . q I . 2 1 ~2 

^n\q) = : sm-Wsm - + sm y„ cos^ - 

/i sm T/„ 2 y 2 2 

with j/„ — (7rn/2/i)(7r — /i) emerge successively sX 
7r/(l + 1/n) from the upper continuum boundaryE3'L3 
All these excitations carry some spectral weight, at 
least for finite N , but most of the spectral weight in 
5„„(g,w) is transferred from the shrinking continuum to 
the top branclh-|the one already present in the WCCF 
reconstruction.c3 

At the endpoint of the Luttinger liquid phase Jz/t = 4, 
the continuum states have been replaced by a series of 
branches e„(q) = {2t/n){\ — cosq), n — 1,2, . . ., all the 
spectral weight is carried by the top branch (n = 1), 
and the dynamic charge structure factor reduces to the 
single-mode form 



S'„„((j,w) = 7r^(5(q)(5(tj) + -6 [u 



T • 2 9 
Jz sm - 
2 



In the framework of the recursion method applied to the 
exact finite-size ground state (|3.2| ) , this simple result fol- 
lows from a spontaneously terminating continued fraction 
with coefficients Ai(g) = sin^(g/2), A2(g) = 0. 

The dynamically relevant charge excitation spectrum 
of i?t-j, which has an even more complex structure, will 
be presented in a separate study. In this case, exact 
results exist onlw at one point {J /t = 2) in the strong 
coupling regime. r 
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E. Spin structure factor 



F. Spin dynamics 



The long-distance asymptotic behavior of the t-J spin 
correlation function in the Luttinger liquid phase was 
predicted t©,bc governed by two leading power-law terms 
of the form&tl 



1 



cos(7rm/2) 

^i)p/4+l 



(3.5) 



where rjp is the charge correlation exponent discussed pre- 
viously. The open circles in Fig. |(a) depict the spin 
structure factor Szz{q)t-j for J/t = 0+ of a system with 
A'^ = 56 sites as inferred via numerical Fourier trans- 
form from the results for the spin correlation function 
presented in Sec. II. The two asymptotic terms of (3.5) 
are reflected, respectively, in the linear behavior at small 
q and in the pointed maximum at q = 7r/2. The lat- 
ter turns into a square-root cusp as N 00. The ex- 
trapolated maximum is Szz{'^/'2)t-,j = 0.28(1) (indicated 
by a -I- symbol). The extrapolated slope at = is 
Szz{q)t-j/q = 0.0847(20). The observed smooth mini- 
mum at q = TT suggests that Szz{q)t-j, unlike S'„„(g)t_j, 
has no singularity there. The extrapolated value is 
Sz.(Tr)t-j - 0.127019(2). 



The predictions of (3.5) that the linear behavior in 
Szz{q)t-j at small q persists throughout the Luttinger 
liquid phase and that the cusp singularity at q = 7r/2 
weakens with increasing J/t and disappears at the on- 
set of phase separation are consistent with our result for 
J/t = 3.2, plotted in Fig. ^(b). The open circles suggest 
a smooth curve which rises linearly from zero at q — Q. 
The smooth extremum at g = tt has turned from a min- 
imum at J/t = 0^ into a maximum at J/t = 3.2. 

The solid line in Fig. |^(a) represents Szz{q)t-.h for the 
free-fermion case Jz/t = 0"*" as obtained from Fourier 
transforming (2.9). It differs from the corresponding 
t-J result (o) mainly in three aspects: (i) the rise 
from zero at small q is quadratic instead of linear, re- 
flecting non-singular behavior at q = 0, i.e. the ab- 
sence of a non-oscillatory power-law asymptotic term in 
{S[S[_^^)t-j/, (ii) the singularity at q = 7r/2 is divergent: 
~ |g — 7r/2|^^/^; (iii) the smooth local minimum at q = tt 
has a slightly higher value, Szz{'^)t-.j^ — 0.129. 

Over the range of the Luttinger liquid phase, the 
asymptotic term in {S[ Si^^)t-j^ which governs the 
singularity in Szz{q)t-j^ at q ~ 7r/2 is of the form 
~ B2 cosiTTm /2) /irfp/^ . As in the t-J case, the singular- 
ity weakens gradually and then disappears at the transi- 
tion point, Jz/t = 4. The finite- A^ result of Szz{q)t-j^ at 
Jz/t = 4, (•) in Fig. ^(b), indeed suggests a curve with 
no singularities. This is confirmed by the exact result. 



Szz{q) 



cos q) , 



(3.6) 



inferred from the exact ground-state wave function ( |3.2| ) 
for N 00. It reflects a spin correlation function which 
vanishes for all distances beyond nearest neighbors. 



Under mild assumptions, which have been tested for 
Ht-j^ at Jz/t = 0+, the following properties of the dy- 
namic spin structure factors Szziq,^) of Ht-j or Ht-j^ 
can be inferred from the singularity structure of Szz{q)'- 
(i) The excitation spectrum in Szz{q,^) is gapless at 
q ~ tt/2. (ii) The spectral-weight distribution at the 
critical wave number q — tt/2 has a singularity of the 
form: 



Sz 



V2' Jt- 



In the weak-coupling limit (rjp = 2), this yields ^ uj '^/^ 
for Ht-j- and ~ w"^/^ for Ht-j- In both cases, the in- 
frared exponent increases with increasing coupling. A 
landmark change in ^^^(Tr, w) occurs at the point where 
the infrared exponent switches sign (from negative to 
positive). In the t-Jz case this happens for rjp — 8 and in 
the t-J case for rjp = 4. According to the data displayed 
in Fig. H this corresponds to the coupling strengths 
Jz/t — 3.6955 . . . and J/t ~ 2.3, respectively. 

The dynamic spin structure factor Szz{q,'-^)t-j^ as ob- 
tained via the recursion method combined w ith a strong- 
couphng continued-fraction (SCCF) analysi£3£a is plot- 
ted in Fig. I as a continuous function of lu and a discrete 
function of g = 2TTm/N, m = 0, . . . , A^/2 with A^ = 12 for 
coupling strengths Jz/t — 0+,2,3,4. This function has 
a non-generic {q it — q) symmetry, which obtains for 
the dynamically relevant excitation spectrum and for the 
line shapes, but not for the integrated intensity. In the 
weak-coupling limit, Jz/t = 0"*", the spectral weight in 
Szzily'-^) is dominated by fairly well defined excitations 
at all wave numbers. The dynamically relevant disper- 
sion is I cos g I -like. 

With Jz/t increasing toward the endpoint of the Lut- 
tinger liquid phase, the following changes can be observed 
in Szz{q,'-^)- The peaks at q ^ 7r/2 gradually grow in 
width and move toward lower frequencies. The |cosg|- 
like dispersion of the peak positions stays largely intact, 
but the amplitude shrinks steadily. The central peak at 
the critical wave number q ~ -k /2 starts out with large 
intensity and slowly weakens with increasing coupling. 
Between Jz/t = 3 and Jz/t = 4, it turns rather quickly 
into a broad peak, signaling the expected change in sign 
of the infrared exponent. 

The dynamically relevant dispersion of the dominant 
spin fluctuations as determined by the peak positions in 
our SCCF data for Szz{q, w) is shown in Fig. |6|for several 
values of Jz/t. The linear initial rise from zero at q — 7r/2 
is typical of a Luttinger liquid. The amplitude of the 
I cos q I -like dispersion decreases with increasing Jz/t and 
approaches zero at the transition to phase separation. 
At the same time, the line shapes of Szz{q, ^)t-j^ tend to 
broaden considerably. These trends are not shared with 
the t-J spin excitations as we shall see. 

The SCCF analysis indicates that the Luttinger liquid 
phase of the t-J model can be divided into two regimes 
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with distinct spin dynamical properties. For coupling 
strengths < J/t < 1, the function Szz{q,i^)t-j, which 
is plotted in Fig. 0, exhibits some similarities with the 
corresponding t-J^ results. The main commonality is a 
well-defined spin mode at not too small wave numbers 
with a I cos g|-like dispersion. This dispersion is displayed 
in the main plot of Fig. ^ for different J/t- values within 
this first regime of the Luttinger liquid phase. 

However, even in the common features, the differences 
cannot be overlooked: (i) The {q tt — q) symmetry in 
the line shapes of 5'^z(9, w)t-j^ is absent in Szz{q,oj)t-j. 
(ii) The amplitude of the |cosg|-like dispersion grows 
with increasing J/t, contrary to the trend observed in 
Fig. H for the corresponding t-Jz spin dispersion, (iii) The 
gradual upward shift of the peak position in 6*22 (tt, u!)t-j 
is accompanied by a significant increase in line width 
(see inset to Fig. ||). Over the range < J/t < 1.25, 
the trend of the q ~ ir spin mode is opposite to what 
one expects under the influence of an antifcrromagnetic 
exchange interaction of increasing strength, (iv) The in- 
tensity of the central peak in SzziT^/2,i^)t-j is consider- 
ably weaker than in in S'z2(7r/2, w)f_j^ . The peak turns 
shallow and disappears quickly with increasing coupling 
(see Fig. ^ main plot). This observation is in accord with 
the proposed dependences of the infrared exponents on 
the coupling constants, (v) The linear dispersion of the 
dynamically relevant spin excitations have markedly dif- 
ferent slopes above and below the critical wave number 
q = n/2 (Fig. ^, main plot). At long wavelengths the 
spectral weight in Szz{q,^)t-j is conceatrated at much 
lower frequencies than in Szz{q,i^)t-j^S^ 

As the coupling strength increases past the value J/t ~ 
0.75, the spin modes which dominate Szziq,^)t-j in the 
first regime of the Luttinger liquid phase broaden rapidly 
and lose their distinctiveness. There is a crossover region 
between the first and second regime, which roughly com- 
prises the coupling range 1 < J/t ^ 2. Over that range, 
the spin dynamic structure factor tends to be governed 
by complicated structures with rapidly moving peaks. 

At the end of the crossover region, a new type of spin 
mode with an entirely different kind of dispersion has 
gained prominence in Szziq,oj)t-j, and it stays dominant 
throughout the remainder of the Luttinger liquid phase. 
This is illustrated in Fig. |l^ for three J/t- values in the 
second regime of the Luttinger liquid phase. The disper- 
sion of these new spin modes gradually evolves with in- 
creasing coupling strength as shown in the inset to Fig. ||. 
Note that the frequency has been rescaled by J both here 
and in Fig. At J/t < 2.0 the dispersion has a smooth 
maximum at g = tt and seems to approach zero linearly 
as g ^ 0. As J/t increases toward the transition point, 
the peak positions in Szziq,'^)t-j gradually shift to lower 
values of w/J, most rapidly at q near tt. 



IV. PHASE SEPARATION 

The transition from the Luttinger liquid phase to a 
phase-separated state in iJt_j^ takes place at J^/t = 4. 
The equivalent XXZ model undergoes a discontinuous 
transition to a state with ferromagnetic long-range or- 
der at the corresponding parameter value (J||/J± = 1). 
The ground state at the transition is non-critical and de- 
generate even for finite N. The XXZ order parameter, 
M = iV-i J2i Sf, commutes with Hxxz- 

Notwithstanding the exact mapping, the transition of 
Ht-j^ at J^/t = 4 is of a different kind. Only one of 
the -|- 1 vectors which make up the degenerate XXZ 
ground state at J\\/J± = 1 is contained in the invari- 
ant subspace that also includes the t-Jz ground state. 
The other vectors correspond to t-Jz states with differ- 
ent numbers Ne of electrons. The t-Jz ground state at 
Jz/t = 4 for fixed = N/2 is non-degenerate and rep- 
resented by the wave function \(f)o) as given in (^.2|). 

The fully phase-separated state as represented by the 
wave function 

1 ^ 

\^,)^-=J2\h,h + l,...,h+N/2-l) 

x{|tiT...)±UTi...)} (4.1) 

has an energy expectation value at Jz/t ~ 4, {Ei) = 
—t{N — 2), which exceeds the finite- A^ ground-state en- 
ergy, Eq = —tN, pertaining to |^o)- However, by com- 
paring the Jz-dependence of the energy expectation val- 
ues (per site) of the two wave functions |0o) a-nd \(t>i), 

in the vicinity of the transition, Jz/t — 4(l+e), we obtain 

_ N^oo e 
eo-er -, 

which implies that the two levels cross at Jz/t = 4 in the 
infinite system. 

The transition to phase separation in iJf.j^ is charac- 
terized by the charge and spin order parameters, 

1=1 1=1 

Neither operator commutes with iJt_ j, . The phase- 
separated state of Ht_j^ is characterized, for A^ 00, 
by a broken translational symmetry, (Qp) 7^ 0, and a 
broken spin-fiip symmetry, (Qa-) 0. 

In the t-J model, the transition to the phase-separated 
state, which takes place at J/t ~ 3.2, produces charge 
long-range order, (Qp) ^ 0, but is not accompanied by 
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the onset of spin long-range order, {Qa) = 0. The sim- 
ilarities in the charge correlations and the differences in 
the spin correlations of the two models are evident in the 
finite-size static charge and spin structure factors. 



A. Charge structure factor 

The vanishing charge correlations in the finite-size t-J^ 
ground state at the onset of phase separation {Jz/t = 4) 
is reflected in the flat charge structure factor Snnio) as 
shown in Fig. 11(a). The corresponding t-J result for 
J /t ~ 3.2 as shown in Fig. 11(b) indicates that correlated 
charge fluctuations do exist at the transition. 

With the exchange coupling increasing beyond the 
transition point, the charge structure factors of the two 
models become more and more alike and reflect the char- 
acteristic signature of phase separation. Phase separa- 
tion is associated with an enhancement of Snn{q) in the 
long- wavelength limit. Because of charge conservation, 
this enhancement is manifest, in a finite system, not at 
q = but at q = 2'k/N. It is conspicuously present in 
the data for couplings J^/t = 4.5 and J/t — 3.5, not far 
beyond the transition point. 

The charge correlation function for the fully phase sep- 
arated state, as represented by the wave function (4.1), is 
a triangular functionjlj {nini+m) 



1/2- \m\/N, |m| < 
N /2. This translates into a charge structure factor of the 
form 



Snn{q) 



N 



4 ^9^0 ■ Ar(l-cosg) 

as shown (for N ~ 12) by the full diamonds in Fig. |l^. 
This function vanishes for all wave numbers q — 2t:1/N 
with even I and increases monotonically with decreasing 
odd The data in Fig. |ll| suggest that the phase sep- 
aration is nearly complete before the exchange coupling 
has reached twice the value at the transition. In the t-Jz 
case, we already know that complete phase separation is 
established (for N oo) right at the transition. 

B. Spin structure factor 

The extremely short-ranged spin correlations in the 
t-Jz ground state (3.2) for N oo are reflected by the 
static spin structure factor (3.6). For finite N the spin 
correlations at distances |n| > 2 do not vanish identi- 
cally. An exponential decay is observed instead with a 
correlation length that disappears as — + oo. Hence the 



1 + cos{Nq/2) 



(1 - <5,,o), (4.2) 



difference between (3.6) and the finite- A data depicted 
in Fig. 12(a) (•). The t-J spin structure factor near the 
transition (J/t ~ 3.2) has a similar q-dependence ex- 
cept at small q, where it tends to zero linearly instead of 
quadratically. 

Whereas the charge structure factors of the two models 
become more and more alike as the exchange coupling in- 
creases in the phase-separated state (Fig. ni|), divergent 



trends are observed in the respective spin structure fac- 
tors, on account of the fact that the t-Jz model supports 
spin long-range order, and the t-J model does not. 

The fully phase-separated state of the t-Jz model is 
at the same time fully Neel ordered. The spin cor- 
relation function in the state (|4.l[) reads {S[S[_^j-^) — 
i(-l)™(l/2-|m|/A), |m| < A/2 and the corresponding 
spin structure factor has the form 

o / ^ ^ <: 1 - cos[A(7r - g)/2] ^ , , ^ ^. 
SM = -,S,..+ 4A[l-cos(.-,)] (^-'^-)- (4-3) 



The function (4.3) vanishes (for even A/ 2) a t all wave 
numbers q = 2ttI/N with even Z, just as (4^) did. The 
exception is the wave number q — t:, where Szz{q) as- 
sumes its largest value. 

The t-J spin structure factor evolves quite differently 
in the presence of increasing phase separation as is illus- 
trated in Fig. 12(b). The electron clustering produces in 
this case the Heisenberg antiferromagnet, whose ground 
state is known to stay critical with respect to spin fluctu- 
ations. The spin structure factor of that model is known 
to be a monotonically increasing function of q, which 
grows linearly from zero at smalLp and (for A —^ oo) 
diverges logarithmically at q = tt.ES 



C. Spin dynamics {t-J model) 

The charge long-range order in the phase-separated 
state freezes out the charge fluctuations in both mod- 
els, and the accompanying spin long-range order in the 
t-Jz model freezes out the spin fluctuations too. What 
remains strong are the spin fluctuations in the t-J model. 

At the transition to phase separation {J/t ~ 3.2), the 
q ~ TT spin mode in Szziq,i^)t-j does not go soft. How- 
ever, the gradual electron clustering tendency in con- 
junction with the continued strengthening of the anti- 
ferromagnetic exchange interaction brings about a soft- 
ening in frequency and an enhancement in intensity of 
the order-parameter fluctuations associated with Nccl or- 
der. Both effects can be observed in the reconstructed 
dynamic spin structure factors at J/t = 3.25,4.0,5.0 as 
shown in Figs. 11(c), 13(a), and 13(b). 

A close-up view of the gradual transformation of the 
q = n mode is shown in Fig. 14(a). For sufficiently strong 
exchange coupling, the function Szzi^^, ^)t-.] will be char- 
acterized by a strongpiiC. nonintegrable infrared diver- 
gence, ~ V— In w/ti;,Eil which characterizes the order- 
parameter fluctuations of the ID s = 1/2 XXX anti- 
ferromagnet. 

Figure 14(b) shows the gradual change in line shape 
and shift in peak position of the function Szz{'^/'2,i^)t-j 
in the phase-separated state. The peak, which starts out 
relatively broad at the transition, shrinks in width, loses 
somewhat in intensity, and moves to a higher frequency. 
For J/t ^ 5.0 it settles at cu/J ~ tt/2 in agreement with 
the lower boundary, uJL{q) — {'^J/'^) \ sing|, at g = 7r/2 of 
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the 2-spinon continuum. The width has shrunk to a value 
consistent with the width of the 2-spinon continuum at 
that wave number. 

In the inset to Fig. 14 we show the evolution of the 
dynamically relevant dispersion for Szz{Q,^)t-j in the 
phase-separated state, as determined by the peak posi- 
tions of our data obtained via SCCF reconstruction. The 
dashed line represents the exact lower threshold of the 2- 
spinon continuum. The shift of the peak positions in our 
data is directed toward that asymptotic position at all 
wave numbers for sufficiently large J/t. 
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FIG. 1. Static charge structure factor at T = of (a) the 
t-Jz model and (b) t-J model in the Luttinger liquid phase. 
Results extracted from the ground-state wave function deter- 
mined numerically for a system of A'^ = 12 sites. 



FIG. 2. Main plot: Renormalized bandwidth of the dy- 
namically relevant charge excitations in the weak-coupling 
regime of the Luttinger liquid phase of the t-J and t-Jz mod- 
els. Inset: Charge velocity in the two models over the full 
range of the Luttinger liquid phase. The open symbols rep- 
resent weak-coupling continued- frac tion data and the solid 
lines represent the exact expression (3.3). The full circles are 
finite-chain data from Ref. n. 



FIG. 8. Dynamically relevant dispersions of the excitations 
dominating the dynamic spin structure factor Szz{<l,i^) at 
T = for t = 1 and difi'erent values of J in the first regime 
(main plot) and the second regime (inset) of the Luttinger liq- 
uid phase of the t-J model. The symbols, which are smoothly 
interpolated by solid lines, represent the peak position of re- 
sults such as shown in Figs. M and nol 



FIG. 9. Line shape at g = tt (inset) and q — ■k/2 (main 
plot) of the dynamic spin structure factor Szz{q,^) at T — 
for 4=1 and various values of J in the first regime of the Lut- 
tinger liquid phase of the t-J model. The results are obtained 
by the same method as those in Fig. ^ 



FIG. 3. Inset: Infrared exponent j3p as defined by (3.4) 
in the weak-coupling regime of the Luttinger liquid phase 
of the t-J and t-Jz models. Main Plot: Inverse square of 
the charge correlation exponent for both models over the full 
range of the Luttinger liquid phase. The open symbols rep- 
resent weak-coupling continued-fraction data, the solid lines 
represent the exact expression (3.1), and the short-dashed line 
the same expression with J/3.2t substituted for Jz/4:t. The 
full circles are the finite-chain data from Ref. ^ 



FIG. 4. Static spin structure factor at T = of the t-Jz 
and t-J models (a) in the free-electron limit and (b) at the 
transition to phase separation. The data for Jz — 0' ^ ar e 
calculated via numerical Fourier transform of expressio n (|2.9|) . 
The data for J = 0^ are derived from expression ( |2.7[ ) as 
explained in the text. The remaining results are extracted 
from the ground-state wave function determined numerically 
for systems of A'' = 12 sites. 



FIG. 5. Dynamic spin structure factor Szz{q,iL>) at T = 
in the Luttinger liquid phase of the t-Jz model. The results 
for t — 1 and four different values of Jz axe obtained via 
strong-coupling continued-fraction reconstruction from the 
coefficients Ai,...,A6 and an unbounded gap terminator 
(Refs. |^j28|). The Afc's are extracted from the ground-state 
wave function for a system of = 12 sites. 



FIG. 6. Dynamically relevant dispersions of the excitations 
dominating the dynamic spin structure factor Szz(g,tj) at 
T = for t = 1 and different values of Jz within the Lut- 
tinger liquid phase of the t-Jz model. The symbols, which 
are smoothly interpolated by solid lines, represent the peak 
positions of results such as shown in Fig. ^ 



FIG. 10. Dynamic spin structure factor Szziq,(j->) at T = 
in the second regime of the Luttinger liquid phase of the t-J 
model. The results for t = 1 and three different values of J 
are obtained by the same method as those in Fig. H. 



FIG. 11. Static charge structure factor at T = of (a) 
the t-Jz model and (b) the t-J model in the phase-separated 
state. Results extracted from the ground-state wave function 
determined numerically for systems of A*' = 12 sites. 



FIG. 12. Static spin structure factor at T = of (a) the 
t-Jz model and (b) the t-J model in the phase-separated state. 
Results extracted from the ground-state wave function deter- 
mined numerically for systems of A' = 12 sites. 



FIG. 13. Dynamic spin structure factor Szz{q,(^) at T = 
in the phase-separated state of the t-J model. The results for 
t = 1 and two values of J are obtained by the same method 
as those in Fig. |[ 



FIG. 14. Line shape of the dynamic spin structure fac- 
tor (a) ^^^(Tr,^;) and (b) Szz{t^/2,u)) of the t-J model in 
the phase-separated state. The results for t — 1 and vari- 
ous values of J are obtained by the same method as those 
in Fig. ^ Inset: Dynamically relevant dispersions of the 
excitations dominating the dynamic spin structure factor 
Szz{q,(^) at T = for t = 1 and difi'erent values of J in the 
phase-separated state of the t-J model. The symbols which 
are smoothly interpolated by solid lines, represent the peak 
position of results such as shown in Fig. n3l 



FIG. 7. Dynamic spin structure factor Szz{q,(^) at T = 
for t = 1 and two values of J in the first regime of the Lut- 
tinger liquid phase of the t-J model. The results are obtained 
by the same method as those of Fig. 5. 
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